Anticipated synchronization: a metaphorical linear view 
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We study the regime of anticipated synchronization recently described on a number of dynamical 
systems including chaotic and noisy ones. We use simple linear caricatures to show the minimal 
setups able to reproduce the basic facts described. 
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I. INTRODUCTION 

The coupling of dynamical systems can lead to the syn- 
chronization of their outputs [lj . Recently, attention has 
been given to the a peculiar phenomenon introduced by 
Voss 0, 13- 0, El where one system synchronizes not with 
the present state but with some future state of another 
system. This regime, called anticipated synchronization, 
has been demonstrated theoretically and numerically in 
disparate dynamical systems including discrete or con- 
tinuous, chaotic and noisy excitable, and quasi-periodic 
ones |E S IE H El • Experimental results have consid- 
ered either electronic circuit implementations of the dy- 
namical equations 5, 9] or lasers running in the chaotic 
regime both in one-directional an d bidirectional cou- 

pimg in eiii m. 

From the viewpoint of an external observer, the dy- 
namics of anticipated synchronization can be seen as 
if one system is forecasting the state of the other. It 
is unclear, however, under which general conditions a 
given dynamical system would or not exhibit such dy- 
namics. The phenomenon in itself is rather counterin- 
tuitive because, depending on the setting, mixes notions 
of synchrony and order with the dynamics of chaos and 
disorder. Furthermore, it is not immediately apparent 
whether or not nonlinearities are essential for the pro- 
cess. This paper is dedicated to answer these concerns, 
by analyzing the minimal setup able to exhibit the more 
fundamental aspects of this phenomena. 

The paper offers three perspectives covering from the 
simplest scenario to the more complex ones. Section II 
is a brief overview of anticipated synchronization as de- 
scribed in recent work. Section III discusses a toy model 
describing the trivial situation of two "particles" moving 
with uniform trajectory in which one attempts to closely 
follow the other using the scheme described in the an- 
ticipated synchronization literature. These particles can 
be seen as an special case of the dynamical system pre- 
sented in Section IV, where the case of two coupled maps 
is analyzed from a dynamical systems perspective and the 
condition for anticipated synchronization derived analyt- 
ically. The same system is revisited in Section V, but 
now from a Control Systems point of view, both for con- 
tinuous and discrete cases. The paper closes by listing 



the most relevant conditions one would expect to see in 
all cases of anticipated synchronization. 



II. OVERVIEW OF ANTICIPATED 
SYNCHRONIZATION. 

Two different schemes have been proposed in order to 
achieve anticipated synchronization @|. Both schemes 
use, in a way or another, delay lines which allow forecast- 
ing of a master trajectory by a siave identical system. 

The first scheme uses the technique of complete re- 
placement. It considers a dynamical system, x(t) (known 
as the master system) whose dynamics involves a delayed 
feedback of the form: 



x(t) = -ax(i)+f(x(t-T)), 



(1) 



being a > a constant, and f a given, generally non- 
linear, function. The dynamics of the siave system, y(t), 
is obtained by a similar equation in which the delay 
term has been completely replaced by the master sys- 
tem. Namely, the evolution equation for y(t) is: 



y(t) = -oy(t)+f(x(t)) 



(2) 



It is easy to see that the manifold y(t) = x(t + r), in 
which the slave anticipates by a time r the actual output 
of the master, is a stable solution of the dynamical equa- 
tions This follows readily from the (exact) evolu- 
tion equation A(t) = — aA(i) for the delayed difference 
A(i) = y(i) — x(t+r). This result is independent of both 
the function f(x) under consideration, or the arbitrarily 
large delay time r. The structural stability of the asymp- 
totic solution has been demonstrated by implementation 
in electronic circuits 0. Of course, the result is more 
remarkable when the dynamics of the master has a high 
degree of unpredictability, i.e. by choosing f(x) and r 
such that the dynamics of the master is chaotic. In this 
scheme, and beyond the mathematical result, the antici- 
pation mechanism can be understood from the fact that 
the dynamics of the master at time t influences its own 
dynamics at a time r later, whereas it enters the dynam- 
ics of the slave immediately at time t. In other words, it 
is as if the slave system y(t) has anticipated knowledge 
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of an essential part of the dynamics of the master system 
x(t). 

The second scheme is one that includes only delay in 
the slave dynamics |2|. This is defined by the following 
dynamical equations for the master and slave systems: 



*(*) = f(x(t)) 

m = f(y(*)) 



•K[x(t)-y(t-r)], 



(3) 



where f (x(i)) is an arbitrary function and K is a coupling 
strength matrix. It is easy to show that the anticipated 
synchronization manifold y(t) = x(i + r) is again a solu- 
tion of this set of equations. At variance with the case of 
complete replacement, its stability can only be fulfilled in 
a limited range of parameters for r and K. Despite this 
restriction, however, we believe that this method of an- 
ticipated synchronization deserves more attention than 
that of complete replacement because the anticipation 
time t is included as an external parameter and does not 
influence the dynamics of the master system. Therefore, 
in principle, any dynamical system can be predicted by 
using this scheme. Again, the anticipation is more rele- 
vant when the dynamics of the master is unpredictable 
by other simple means. 

Recently an extension of the scheme was introduced to 
include anticipation in non-autonomous systems |9j, [10( . 
Specifically, it has been considered the following set of 
equations: 



x(t) = f(x(t))+I(t) 

y(i) = f(y(i))+I(t)+K[x(t)-y(i-r) 



(4) 



where the new term l(t) represents an external input 
acting on both, master and slave, systems. Although 
the anticipated manifold y(t) — x(t + r) is no longer an 
exact solution of the previous equations (except in the 
case of a periodic forcing I(f + r) = I(i)), it has been 
shown that several features of the dynamics of the master 
can indeed be predicted by the slave. For instance, it is 
possible to predict the peaks fired by an excitable system 
subjected to a random external forcing [9|,ll0|. We stress 
that in this case the random forcing induces peak firing 
at uncorrelated and unpredictable times. 

In this second scheme, Eqs. Il.'il II) . the actual mech- 
anism leading to synchronization is much more elusive 
and, despite the wide variety of work, it is still unclear 
which are the relevant conditions and requirements for 
two dynamical systems to exhibit this type of anticipated 
synchronization. For instance, one would like to under- 
stand whether or not nonlinear aspects of the dynamics 
are needed for the systems to exhibit anticipated synchro- 
nization. In the same line to what extent anticipation can 
be arbitrarily long is of relevance for practical purposes. 



III. TWO PARTICLES FOLLOWING EACH 
OTHER 

The intention in this section is to undress the antici- 
pated synchronization models of irrelevant aspects, to be 



able to identify the essential mechanisms at play. The 
goal is to have a system homologous to the more general 
ones described in the anticipated synchronization litera- 
ture. We first look at the simplest linear one, thus we 
choose a dynamical system of two particles moving uni- 
formly in a one-dimensional space. We consider a particle 
following a uniform motion 



x(t) = v 



(5) 



where x(t) is the position of the particle at time t and v 
its velocity. Following the scheme in anticipated synchro- 
nization to forecast (and "synchronize" with) the position 
of the master particle we consider another ( "slave" ) par- 
ticle whose trajectory is given by 



y{t) = v + K{x{t)-y{t-T)) 



(6) 



where K is the strength of the coupling between the 
master particle (Eq. (JSJ) and the slave (Eq. ||BJ). To 
achieve the anticipation, the solution y(t) — x(t + r) 
has to be a stable solution. The stability can be studied 
analyzing the behavior of the delayed difference A(i) = 
y(t) — x(t + t). We obtain 



A{t) = -KA(t - t) 



(7) 



The condition for local stability of the solution A(t) = 
(y(t) = x(t + r)) is given by pi lil) 



Kt<* 
2 



(8) 



This indicates that in order to observe an anticipated so- 
lution of the slave particle within this scheme, the prod- 
uct of the coupling with the anticipation time r has to 
be smaller than a certain value. Larger values produce 
an over-correction and the slave is unable to anticipate. 
This expression qualitatively reproduces the numerical 
results obtained by several researchers. 

It is interesting to note that if the two systems are not 
identical, e.g., the velocities are different, then 



+ K(x(t)-y(t-r)) 



(9) 



then we obtain the same stability condition. How- 
ever the anticipation is not perfect in the sense 
that the slave anticipates with a constant mismatch 

v — v 

y(t) = x(t + t) H — — . The mismatch decreases with 

K 

increasing coupling constant K. 

More generally one could obtain a stability equation of 
the from 



A(t) = AA(t) - KA(t 
leading to the conditions 0, 0] 

K < \A\ 



r) , 



or, 



K > \A\ and r < 



cos-\A/K) 



(10) 



(11) 



(12) 



(K 2 - A 2 ) 1 / 2 ' 
where the principal value is taken (0 < cos~ 1 (A/K) < 

7T). 
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IV. A DYNAMICAL SYSTEMS PERSPECTIVE: 
TWO COUPLED MAPS 

Let us consider the same two particles of the previous 
section but now moving with the more general dynamics 
given by the following coupled maps: 



x k+1 = ax k + a 



mod(m) 



Uk+i = ay k + a + j(x k - y k - n ) mod(ro), 



(13) 



where a > and a are constants. If a > 1 the x k 
form a chaotic map, whereas a = 1 leads to a quasi- 
periodic map and for a < 1 the map converges to 
x k = a/(l — a). The necessary n + 2 initial conditions 
are the set (x ,y_ n ,y_ n+1 , ■ ■ ■ ,y-i,y ). 

We are interested in the possibility of the coupled maps 
leading to anticipated synchronization, i.e. in having 



Vk = Xn+k 



(14) 



as an asymptotic solution for the maps. To this end, 
we consider the map satisfied by the delayed difference 

fa = x n+k - Vk'- 



Sk+i = a8 k - 7<5fc_ n . 



(15) 



An analysis complementary to ours can be found in 



The ansatz S k = X k leads to 



n+l 
z=l 



(16) 



The (complex) constants Ci, . . . , C n +i are determined by 
the initial conditions. The A^, i = 1,. .. ,n + 1 are the 
solutions of the polynomial equation: 



A n (A~a) + 7 = 



(17) 



A necessary and sufficient condition for anticipated syn- 
chronization to hold and to be asymptotically stable 
is that lim 5 k = 0, or equivalently |A;| < 1, Vi = 

k — >oo 

1, . . . , n + 1. We analyze now the range of validity of 
this condition as a function of the parameters a, n and 
7. It turns out that anticipated synchronization occurs 
if 7 G (a — 1, 7 C ), where the value of 7 C depends on a and 
n. We consider three cases: 

1.- a < 1. We treat for completeness this case, al- 
though it is the less interesting from the practical point 
of view, since the maps tend asymptotically to a fix point 
(in practise, after a finite number of steps). Although a 
full analytical solution does not seem to be available in 
this case, it is possible to obtain the asymptotic behavior 
in the case of large anticipation n: 



7c ^l-a + — , a<l, 



(18) 



Notice that anticipation is possible for arbitrarily large 



2.- a = 1. In this case, the dynamics is quasi-periodic. 
The interval of anticipated synchronization can be found 
analytically as: 



7c 



2 sin 



2(2n+T 



(19) 



7?- , the same as the continuous 

2n ' 



In the limit n — > 00, 7 C 

case considered in the previous section. Again, anticipa 
tion is possible for arbitrarily large n. 

3.- a > 1. This is the more interesting case since the 
maps are chaotic. It turns out that the condition |A;| < 1 
is only satisfied for a < 1 + i. Alternatively, for fixed 
a > 1 the maximum anticipation time is 



1 



1 



a > 1 



(20) 



where [x] denotes the integer part of x. Anticipated syn- 
chronization is found for 7 E (a — 1,7c) but only for 
n < n max . In this case of a > 1 we have not been 
able to find any analytical approximation or asymptotic 
expression and the values of 7 C need to be computed nu- 
merically. The previous results are summarized in Fig^ 

Let us now analyze the predictability of the map, tak- 
ing a = 1.1, n = 6, for which the synchronization limits 
are 7 £ (0.1,0.19395). First, for visualization purposes, 
we transform the output of the master map into a series 
of random spikes by defining a new "firing" sequence u k 
as u k = 1 if X}~ > u and u k = if x k < u, with u a 
convenient threshold value. Similarly for the slave map, 
we define v k — 1 if y k > u and v k = if y k < u. In the 
case of a — 1.1 the timing of the pulses does not have 
any regularity. As shown in Fig. [21 the power spectrum of 
the u k signal is flat, showing the absence of any preferred 
time scale, the series is highly unpredictable. 

In Fig. |31 we show that, in accordance with the previ- 
ous analysis, each pulse of the master system u k is an- 
ticipated n = 6 units of time by a pulse in the slave 
system v k - n for a coupling 7 = 0.16 G (0.1,0.19395) 
(middle panel). For too large or too small values of the 
coupling 7, the synchronization is lost. The existence of 
a minimum and maximum value for the coupling in or- 
der to have anticipated synchronization, which appears 
here as a property of simple linear maps, has also been 
observed in more complex chaotic and excitable non- 
autonomous systems [l(|. These maps can be seen as 
well as a case of integrate and fire dynamics, often used 
as caricatures of neural systems. 



V. A CONTROL SYSTEMS PERSPECTIVE 

From an engineer viewpoint, the slave dynamical sys- 
tem following a master, shown in the previous section, is 
seen as a cascade of control system blocks. In this section 
anticipated synchronization is analyzed from this Control 
System perspective. In both, the continuous and discrete 
time systems, x(t) will be seen as the output of a first 
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FIG. 1: The anticipated synchronization regime is achieved 
for values of 7 and n in between the upper (solid) line and the 
lower (dashed) line. For a — 1 the solid line is the analytical 
result Eq. H19I . whereas for a = 0.5 and a = 1.1 the dots 
have been computed numerically and the solid line is a guide 
to the eye. 
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FIG. 3: Pulses Uk (solid lines and circles) and Vk- n (dashed 
lines and diamonds), generated by the iteration of Eq.ED with 
a — 1.1, a = 0.1, n — 6, for three values of 7 . Notice that 
anticipated synchronization (coincidence of the Uk and Vk-n 
pulses) occurs only for the intermediate value of 7 (middle 
panel) . 
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FIG. 2: Power spectrum u u (arbitrary units) of the series 
of Uk pulses generated by Eg. 1131 with parameters a — 1.1, 
a — 0.1 and m = 1. 



order open loop system with no feedback, and the slave 
y(t) (controlled signal) represents the output of an iden- 
tical system but with delayed time feedback and driven 
by the master system x(i) (the reference signal). With 
enough loop gain in the feedback loop, the control will 
act as a servo mechanism which minimizes the error be- 
tween the reference (master x(i)) and its own delayed 
output y(i — r). Hence, the internal variable y(t) will be 
a prediction of x(i). 



Two particles following each other seen as a 
servo mechanism 



The example of the two particles following each other 
described by equations (J5J) and @ can be seem as a servo 
mechanism where the position of the controlled particle 
y(t — t) is compared with a reference signal x(t). The 
integrator nature of the loop provides enough gain so 
that in steady state the error will vanish and the two 
signals will match. 

Equation JSJ represents a pure integrator to input v. 
Denoting by X(s) the Laplace transform of x(t) and as- 
suming zero initial conditions for x(t = 0), the Laplace 
transform of x(t) will be sX(s). Also, since the constant 
velocity v is represented as an external step input, its 
Laplace transform V(s) will be equal to v/s. Thus we 
can write 



X(s) 



v 

~2 



(21) 



This represents a transfer function, in the Laplace do- 
main, with a pole at the origin of the s plane. If v is 
constant the output is a ramp x(t) — vt. 

Repeating the same process for Eq. (j^J) we can write 



sY(s) = - + K(X(s) - e- TS Y(s)) , 



from where we obtain 



Y(s) = 



v/s 



Ke 



KX(s) 
+ Ke- T 



(22) 



(23) 



Applying the inverse Laplace Transform to Eq. (|23|) we 
can obtain the temporal evolution of the output for any 
given input. The dynamics of the error in the servo mech- 
anism can also be obtain by subtracting X(s) — Y(s). 

The input/output relation, (transfer function in con- 
trol system terms), between X(s) and Y(s) is given by 
the ratio of these two functions (considering zero external 
input): 



G cl (s) 



Y(s) 
X(s) 



K/t 



(24) 



Equations (|22|) - (|24|l can be represented in a block dia- 
gram, as in Fig. 0] with direct gain G(s) = — and feed- 
back gain H(s) = e~ TS . In more general terms, the closed 
loop transfer function is given by 



G c l(s) 



Y(s) 

X{s) l + G(s)H(s) 



G(s) 



(25) 



It follows that if the loop gain G(s)H(s) S> 1, the 
closed loop transfer function will be given by the feed- 
back term G c i ~ (i?(s)) _1 . For our example this means 
that Y(s) ~ X(s)e TS , and applying the inverse Laplace 
transform y(t) ~ x(t + t). This was the expected result 
that proves that the delayed position of the second par- 
ticle, after some transient, will follow the position of the 




FIG. 4: Block diagram of the continuous-time servo mecha- 
nism system (Eq. 1221 - 1241 ') represented by a transfer func- 
tion in the Laplace domain. 



first one. Likewise, y(t) will be moving ahead of x(t) in 
a predictive manner. 

To perform a stability analysis, the roots of the denom- 
inator s + Ke~ TS = must be found. If we replace s by 
iui a frequency analysis can be obtained. This is usually 
done drawing real and imaginary parts of the transfer 
function with the Nyquist diagram or using the Bode 
plots, representing module and phase separately. The 
term G(iuj) ~ K/ito has a gain whose module decreases 
with co and becomes unity for lo — K, and a constant 
phase equal to — ir/2. The transfer function of the delay 
H(iuj) has a constant module equal to one 



\H{iu)\ = \e- luJT \ = 1 
and a phase that decreases with w 



[CO 



arg|e 



-LOT 



(26) 



(27) 



According to the standard stability criteria for linear 
systems the module of the gain of G(lo)H(ui) must be 
smaller than 1 when the phase crosses the tt boundary 
(meaning the denominator of G c i(s) is zero). Substitut- 
ing ui = K the total phase contribution [tK + tt/2) at 
the unity gain point should be smaller than 7r, recovering 
Eq. © 



t < 



2K 



(28) 



B. Discrete System 

A discrete system similar to the coupled maps de- 
scribed by Eq. (but without the mod(m) constrain) 
can be seen, in engineering terms, as a discrete time con- 
trol systems written as two difference equations: one for 
the master Xk and the other for the slave yk- The equa- 
tions (|13|) without the mod(m) constrain can be rewritten 
for time k 

x k+ i = ax k + a (29) 
y k +i = ay k +a + ~f(x k - y k - n ) (30) 
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FIG. 5: Block diagram of the discrete-time servo mechanism 
system Eq. 1M It represented by a transfer function in the 
Z-transform domain. 

To solve Eqs. (I29H - I|31)[1 . we will use the Z-transform 
[l8j . If Z{xk} is the Z-transform of Xk denoted as X(z) 
and Z{xk-i} is the Z-transform of Xk—i obtained by 
multiplying X(z) by z , applying these rules to both 
equations, we obtain 



X(z) 
Y(z) 



Z{a} 
1 — az~ x 

Z{a} + 1 z- 1 X{z) 
1 - az- 1 - 7z-(™ +1 ) ' 



(31) 
(32) 



where Z-transform of a, Z{a}, depends on the type of 
input a. For our case, a is constant and 



Z(a) 



(33) 



As in the continuous system, we are interested in the 
transfer function between X{z) and Y (z). Any feedback 
system can be viewed as a direct transfer function G(z) 
and a feedback block H(z) whose output will be sub- 
tracted from the reference, to obtain the error signal. In 
our case 



Y(z) 



G 1 {z)G 2 {z)X{z) + G 2 (z)Z(a) 
l + G 1 {z)G 2 (z)H(z) 



(34) 



where Gi(z) = fz~ 1 , G 2 {z) = (l-az" 1 )- 1 , and H (z) = 
z~ n . In the engineering literature, these equations are of- 
ten represented using block diagrams, as shown in Fig. [5] 
The closed loop transfer function in this case is given by 



G c l(z) 



Y{z) 
X(z) 



7_z 



1 



az 



(35) 



Using the same arguments as in the previous subsection, 
if the loop gain Gi(z)G 2 (z)H(z) > 1 we can approxi- 
mate Gd ~ (H(z))- 1 . This yields to Y(z) ~ z n X(z) 



indicating that the output yk predicts Xk by n sampling 
periods. 

To analyze the stability we analyze the conditions 
when the denominator of Eq. I|34l) equals zero. Thus 
to ensure stability the roots of the denominator must lie 
inside the circle given by 

z n (z-a) = -j, (36) 

recovering Eq. (|17fl . 

VI. CONCLUSIONS 

As stated at the outset, the purpose of this paper is 
to undress the models of anticipated synchronization to 
be able to see the essentials of this intriguing dynamics, 
which in summary are: 

i) The phenomenon of anticipated synchronization in 
itself does not rely on nonlinear properties of the dynam- 
ics, indeed, as shown in Sect. IV nonlinearities make the 
anticipation harder or impossible. 

ii) Some aspects of the problem are naturally under- 
stood by looking it from a Control Systems approach, 
where the delayed term are seen as an "error signal" in 
a servo mechanism control loop (Sect. V). 

iii) As shown by the derivation using simple stability 
criteria (Sect. Ill and IV), the boundaries of the antici- 
pated synchronization (i.e., the fundamental diagram in 
Fig. 1 of coupling strength versus the delay) are expected 
to be universal for this kind of systems. 

Overall, these results account for the necessary quali- 
tative features two systems must have in order to exhibit 
anticipated synchronization. By clarifying under which 
conditions the phenomenon would be observable, some 
questions will naturally arise opening new application of 
these ideas. From the arguments discussed here it seems 
that "synchronizing in advance" and "controlling the fu- 
ture state" can be equivalent objectives describing the 
aim of the process of anticipated synchronization. In 
that sense it might be worth to study possible connec- 
tions with systems that learn to master such objectives, 
as for instance the task of a neural net adapting to cap- 
ture a flying object. 
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